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Magneto-transport properties in closed and open loop structures are carefully reviewed within a 
tight-binding formalism. A novel mesoscopic phenomenon where a non-vanishing current is observed 
in a conducting loop upon the application of an Aharonov-Bohm flux (f> and we explore its behavior 
in the aspects of quantum phase coherence, electron-electron correlation and disorder. The essential 
results are analyzed in three different parts. First, we examine the behavior of persistent current 
in different branches of a zigzag carbon nanotube within a Hartree-Fock mean field approach using 
the second quantized formulation. The phase reversal of persistent current in several branches as 
a function of Hubbard interaction is found to exhibit interesting patterns. Our numerical results 
suggest a filling-dependent metal-insulator transition in a zigzag carbon nanotube. Next, we address 
the behavior of persistent current in an ordered-disordered separated nanotube keeping in mind a 
possible experimental realization of shell-doped nanowire which can provide a strange electronic 
behavior rather than uniformly doped nanowires. Finally, we focus our attention on the behavior of 
persistent current in an open loop geometry where we clamp an ordered binary alloy ring between 
two ideal semi-infinite electrodes to make an electrode-ring-electrode bridge. From our investigation 
we propose that under suitable choices of the parameter values the system can act as a p-type or 
an n-type semiconductor. 
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I. INTRODUCTION 


The abundant progress in nano-science and technology 
has stimulated us to fabricate different artificial nano¬ 
structures over the last few decades whose dimensions 
are comparable to and even smaller than the mean free 
paths or wavelengths of electrons such that they can 
move through the samples without randomizing their 
phase memorieJi^. The manifestation of phase coherence 
length is one of the most striking aspects in the meso¬ 
scopic regime which can be obtained by lowering the 
temperature of the samples below sub-Kelvin tempera¬ 
tures since the scatterings due to electron-electron and 
electron-phonon interactions are highly suppressed in this 
temperature regime. Such experimentally accessible low- 
temperatures also make the system energy levels behave 
like discrete states which play a pivotal role in appear¬ 
ing several quantum-mechanical phenomena. Therefore, 
a mesoscopic system which can be modeled as a phase- 
coherent elastic scatterei!^ provides us the opportunity 
to explore various novel quantum-mechanical effects be¬ 
yond the atomic realirP. The experimental realization 
of different quantum-mechanical incidents like, universal 
conductance fluctuations, non-local current-voltage re¬ 
lationship, new Onsager reciprocity relations. Coulomb 
blockade in micro-tunnel junctions^tSl, Anderson local- 
izatiorP, quantum Hall-effecl!^^, Aharonov-Bohm (AB) 
oscillation^iilli^, to name a few, has raised the popular¬ 
ity of the mesoscopic world among the scientists and en¬ 
gineers. Another reason behind this popularity can be 
attributed to the tailor-made geometries which look very 
simple but have high potential from the application per¬ 


spective. In addition, several other fluctuation patterns 
of the conductance are easily reproducible simply by tun¬ 
ing external parameters like, magnetic held, electric held, 
Fermi level, etc. 

The existence of dissipationless current, the so-called 
persistent current, in a mesoscopic normal metal ring 
pierced by an AB hux cf is one of such remarkable ef¬ 
fects which reveals the importance of phase coherence 
of electronic wave functions in low-dimensional quantum 
systems. The alluring question of persistent current in 
normal-metal rings threaded an AB hux was hrst ex¬ 
plored during 1960’^^. Later, in 1983, Biittiker et al. 
have revived the interest regarding this topicli^. They 
have predicted that a small isolated normal metal ring 
penetrated by a slowly varying magnetic hux carries an 
equilibrium current which does not decay and circulates 
within the sample. But its experimental realization was 
somewhat difficult because at that time it was a chal¬ 
lenging task to conhne magnetic hux in such a small re¬ 
gion like nanoscale ring. However, a few years later Levy 
et al^ have given hrst experimental evidence of persis¬ 
tent current in a mesoscopic metallic ring. They have 
observed the oscillations with period (/)o/2 {cfo = chje, 
the elementary hux-quantum) while measuring persis¬ 
tent current in an ensemble of 10^ independent Cu rings. 
Similar oscillations with period <j)o/2 have also been re¬ 
ported not only for an ensemble of disconnected 10^ Ag 
rijig^SIbut also for an array of 10® isolated GaAs-AlGaAs 
rijig^m Later, many other theoreticaP^USSI as well as ex- 
perimentaP®^ attempts have been done to explore the 
actual mechanisms of persistent current in single-channel 
rings and multi-channel cylinders. However, a contro- 
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versy still persists among the experimental observations 
and theoretical estimates of persistent current ampli¬ 
tudes. All the experimentally measured current ampli¬ 
tudes were found to be one and two orders of magnitude 
larger than the theoretically predicted results, except in 
the case of nearly ballistic GaAs-AlGaAs ring^^ where 
00 -periodic persistent currents have been observed with 
amplitude Jq ~ evp/L {vp and L are the Fermi velocity 
and circumference of the ring, respectively), which is very 
close to the value obtained from the free electron the¬ 
ory at absolute zero temperature (T = OK). In a recent 
work, Bluhm et have examined magnetic response 
of 33 individual cold mesoscopic gold rings, considering 
one ring at a time, using a scanning SQUID technique. 
Their results agree well with the theoretically estimated 
valu^^ in a single ballistic rinj^ and an ensemble of 16 
nearly ballistic ringj^. But, the amplitudes of persistent 
current in single-isolated-diffusive gold rings^ are still 
two orders of magnitude larger than the theoretical cal¬ 
culations. This discrepancy initiated intense theoretical 
activity, and it is generally believed that the electron- 
electron correlation plays an important role in the disor¬ 
dered diffusive ringP^H^S An explanation based on the 
perturbative calculation in presence of interaction and 
disorder has been proposed and it seems to give a quan¬ 
titative estimate closer to the experimental results, but 
still it is less than the measured currents by an order 
of magnitude, and the interaction parameter used in the 
theory is not well understood physically. 

The behavior of persistent current in an isolated loop 
geometry enclosing a magnetic flux is highly sensitive to 
the location of Fermi level of the system where the role of 
magnetic flux is essentially to destroy the time-reversal 
symmetry and as a result, the degeneracies among the 
current-carrying states get removed. Depending on the 
Fermi level and the direction of the magnetic flux cur¬ 
rent flows in either direction revealing diamagnetic or 
paramagnetic nature. Since the discovery persistent cur¬ 
rent have been studying in different kinds of mesoscopic 
loop struc tures. Among them single channel mesoscopic 
ring^SSmi were mainly in the focus of attention in spite 
of having topological simplicity they carry a high poten¬ 
tial from the application point of view. Not only single 
isolated rings but also array of such rings have been used 
to study the nature of persistent current. A very few 
works on multi-channel closed loop systems have been 
discussecP^Ell. in case of multi-channel mesoscopic cylin¬ 
ders it has been noticed that the typical single-level cur¬ 
rent gets reduced with increasing number of conducting 
channels M. The correlations in the energy spectrum 
governs that the ratio of the total and single-level cur¬ 
rents is proportional to 

Similar to persistent current in such isolated closed 
loop geometries, the non-decaying ch arge current 
is also observed in open loop systems^BEBj where 
rings/cylinders are coupled to source and drain elec¬ 
trodes. In 1985, Biittiker has introduced a conceptually 
simple and elegant approaclP^ to investigate the effect of 


an electron reservoir on persistent current in a loop pen¬ 
etrated by a magnetic flux. The reservoir is considered 
as a source of dissipation and the inelastic scattering pro¬ 
cesses take place only within this reservoir to which the 
loop is attached via a single current lead while the other 
scatterings occur in the lead are assumed to be elastic. 
So, there exists a complete spatial separation between 
the elastic and inelastic scatterings. The reservoir being 
a sink of electrons maintained the chemical potential of 
the loop to a fixed value. This modifies the statistical 
mechanical treatment of the system to a different ensem¬ 
ble, the grand canonical ensemble where close systems 
(in absence of lead and reservoir) belong to canonical 
ensemble average. This establishes a remarkable differ¬ 
ence between the study of closed and open loop systems 
from the statistical point of view. Experimental verifica¬ 
tion of persistent current in open system was first done 
by Mailly et al. in 199^^. They have examined both 
closed and open systems and detected a periodic signal 
carrying the signature of persistent current which was in 
agreement with the theory. After those pioneering works 
last two decades have witnessed several approaches to 
reveal the properties of persistent currents in open sys¬ 
tems. Unlike the behavior of the other physical quantity 
such as conductance, persistent current in open system 
is sensitive to the direction of the transport current and 
this property is helpful for recognizing this current from 
the other currents (noise) associated with experimental 
measurements. In open systems, persistent current can 
appear even in absence of external magnetic helcP^ as 
the lead-current, the current from the source to d rain, 
plays the role of the driving force. Jayannavar et 
have shown theoretically the flow of persistent current in 
an open metallic loop connecting to two electron reser¬ 
voirs. They have taken the lengths of the two arms of 
the loop unequal which results in a circulating current in 
the loop. Instead of setting the arm-lengths unequal one 
can also introduce any local scatterer anywhere in the 
geometry to establish a circular flow within the loop ge¬ 
ometry in non-equilibrium situation. On the other hand, 
in an equilibrium condition i.e., when the chemical poten¬ 
tials of both leads are same, persistent current can still 
arise due to evanescent modes. In a work by Jayannavar 
et al. the phenomenon of persistent current due to two 
non-classical effects, AB effect and quan tum mechanical 
tunneling has been studied in detaiP^EU, 

Although the studies involving persistent current in 
single-channel rings and multi-channel ring^^ or cylin¬ 
ders, both in the closed and open loop shapes, have al¬ 
ready generated a wealth of literature there is still need 
to look deeper into the problem to address several signifi¬ 
cant issues those have not been well explored before as for 
examples the understanding of the distribution of persis¬ 
tent currents in different channels in presence of electron- 
electron interaction and disorder which inspect the net 
response of the full system and also the sensitivity of 
persistent current on disorder in partially disordered sys¬ 
tems like shell-doped nanowires which can provide a un- 
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familiar electronic behavior rather than uniformly doped 
nanowires. In the present review we essentially concen¬ 
trate on these issues. 

In the first part we address magneto-transport prop¬ 
erties in a zigzag carbon nanotube, formed by rolling up 
a graphite ribbon in the cylindrical fornP^, with its de¬ 
tailed energy band structure in presence electron-electron 
interaction within a Hartree-Fock (HF) mean field ap¬ 
proach using the second quantized formulation. Since the 
isolation of a single layer graphene by Novoselov et alW^ 
intense and diverse research is going on to explore elec¬ 
tronic transport in this system. Graphene, a single layer 
of carbon atoms tightly packed into a two-dimensional 
honey-comb lattice, has drawn attention of scientists in 
various disciplines due to its unconventional and fasci¬ 
nating electronic properties arising particularly from the 
linear dispersion relation around the Dirac points of the 
hexagonal Brillouin zone. These unique properties can be 
understood in terms of the Dirac Hamiltoniarf^ since it 
actually describes the physics of electrons near the Fermi 
level of the undoped material. The carriers in graphene 
effectively behave as massless relativistic particles within 
a low energy range close to Fermi energy and these mass¬ 
less Dirac Fermion^^ evince various phenomena in this 
energy range. The bipartite character of the wonderful 
lattice structure of graphene strongly influences its intrin¬ 
sic properties and makes graphene a wonderful testbed 
for quantum field theory and mathematical physics as 
well as condensed matter theory. In the last few years ex¬ 
tensive studies on persistent current in carbon nanotubes 
have been performed and ma ny i nteresting physical phe¬ 
nomena have been explorecP®^. Persistent current in a 
carbon nanotube is highly sensitive to its radius, chirality, 
deformation, etc. In a recent experiment it has also been 
observed that the Fermi energy of a carbon nanotube can 
be regulated nicely by means of electron or hole doping, 
which can induce a dramatic change in persistent cur- 
renlP^. It is well established that in a conventional multi¬ 
channel mesoscopic cylinder electron transport strongly 
depends on the correlation among different channels as 
well as the shape of Fermi surface. Therefore, we might 
expect some interesting features of persistent current in 
a carbon nanotube due to its unique electronic structure. 
Here we discuss the behavior of persistent current in dif¬ 
ferent branches together with the total current of a zigzag 
carbon nanotube in presenc e of a n AB flux. Most of the 
works available in literatur^^fllinHlZl generally investigate 
magnetic response of the entire system, but a complete 
knowledge of magnetic response in individual branches 
provide much better insight to the problenP^. The phase 
reversal of persistent current in different branches as a 
function of electron-electron correlation is found to ex¬ 
hibit interesting pattern. Our detailed numerical analysis 
suggests a filling-dependent metal-insulator (MI) transi¬ 
tion in a zigzag carbon nanotube. 

In the second part, we concentrate on the behavior per¬ 
sistent current in an ordered-disordered separated nan¬ 
otube considering a possible experimental realization of 


shell-doped nanowires which may provide several un¬ 
usual electronic behavior rather than uniformly doped 
nanowires and have potential applications in nanoscale 
electronic and optoelectronic devices. 

In the last part we investigate persistent current to¬ 
gether with average density of states (ADOS) in an 
open loop geometry where an ordered binary alloy ring 
threaded by a magnetic flux is clamped between two 
ideal semi-infinite metallic electrodes, commonly known 
as source and drain electrodes, followed by the character¬ 
istic properties of an isolated ordered binary alloy ring. 
Inclusion of some foreign atoms in anyone of the two arms 
of the ring provides some interesting patterns in ADOS 
and from our numerical analysis we propose that under 
suitable choices of the parameter values the system can 
act as a p-type or an n-type semiconductor. 

Throughout the review we choose c = e = h = I for 
numerical calculations and restrict ourselves at absolute 
zero temperature. 

II. MAGNETO-TRANSPORT IN A ZIGZAG 
CARBON NANOTUBE 

This section is devoted to reveal the magnetic re¬ 
sponse of interacting electrons in a zigzag carbon nan¬ 
otube enclosing a magnetic flux within a HF mean field 
approach. Following the description of energy spectra for 
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FIG. 1: (Color online). Schematic view of a zigzag graphite 
nano-ribbon with N^. and Ny number of atomic sites along 
the X and y directions, respectively. 

both non-interacting and interacting cases we investigate 
the energy-flux characteristics, persistent current in in¬ 
dividual branches of the system and also the net current 
of the entire system. 

A. The model and the mean field scheme 

We begin by referring to Fig. [^where a graphite nano¬ 
ribbon of zigzag edges is shown. The filled magenta 
(large) and yellow (small) circles correspond to two dif¬ 
ferent sub-lattices, namely, A and B, respectively. 
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and Ny are the number of atomic sites along the x and 
y directions, respectively. In order to elucidate magnetic 
response of a nanotube we roll up the graphite ribbon 
along X direction using periodic boundary condition and 
allow to pass a magnetic flux </> (measured in unit of el¬ 
ementary flux quantum = chje) along the axis of the 
tube as shown in Fig. Our model quantum system 
is illustrated by the nearest-neighbor tight-binding (TB) 
framework which captures most of the essential proper- 





FIG. 2: (Color online). A graphite nanotube threaded by a 
magnetic flux >1>. 

ties of the tube nicelj^^^HZH Xo incorporate the effect of 
electron-electron interaction in the Hamiltonian we em¬ 
ploy a HF mean field approximation. In Wannier basis, 
the Hamiltonian of an zigzag nanotube takes the form, 

H = t 'y ( „ 1,71,0-6 

m,n,<7 

"b ^ln,n,a^m+l,n,a-£■ + t*m,n,o^oi,7i-l-l,o) + h.C. 

m,n 

where, al„_„,, ibln,n,,^gma) is the creation operator for 
an up spin or down spin electron of spin associated with 


in-l,n m+l,n 



FIG. 3: (Color online). Schematic view of different atomic 
sites with their co-ordinates. 

A (B) type of atoms at the position (m,n) and the 
corresponding annihilation operator is denoted by am,n 
{bm,n)- t is the nearest-neighbor hopping integral and U 
is the strength of on-site Hubbard interaction. The co¬ 
ordinates of the lattice sites are denoted by two integers. 


m and n. The site indexing is schematically shown in 
Fig. I for better viewing. The factor 0 {= 2'iT(j)/Nx), the 
so-called Peierl ’pi phase factor, is introduced into the 
above Hamiltonian to incorporate the effect of magnetic 
flux applied along the axis of the tube. 

Decoupling of interacting Hamiltonian: Using the gener¬ 
alized HF approach, the so-called mean-field approxima¬ 
tion, we decouple the TB Hamiltonian into three different 
parts two of which correspond to two different values of 
a, t and with modified site-energies. The decoupled 
Hamiltonian is expressed as, 

J/mf = H^ + H^ + Ho (2) 

where, 

m,n 

m,n 

+ (^^77i-H.rt.t)^m-l-1.77.i) 

m^n 

m,n 

+ ^L,n,l^-rn,n+l,l + h-C.) , 

Ho = 

m,n 

~b (^m-t-l,n,^)) * (^) 

where, and represent the up-spin and down-spin 
Hamiltonians, respectively, ffg is a constant term which 
gives the energy shift. Here, and are the 

number operators associated with the A and B sites, re¬ 
spectively. 

Self-consistent procedure: In order to get the energy 
eigenvalues of the interacting Hamiltonian we go through 
a self-consistent procedure considering initial guess val¬ 
ues of and {n^rn,n,(y)- With these initial values, 

the up and down spin Hamiltonians are diagonalized nu¬ 
merically and a new set of values of (n((j „ and 
are calculated. These steps are repeated until a self- 
consistent solution is achieved. 

Finding the ground state energy: After getting the self- 
consistent solution we determine the ground state energy 
{Eo) at absolute zero temperature (T = OAT) for a par¬ 
ticular filling by taking the sum of individual states upto 
the Fermi level {Ep) for both up and down spin electrons. 
The expression for ground state energy reads, 

Eo = Ei^A; -b Ei^^ -b Ho (4) 

i i 

where, i runs over the states up to the Fermi level. Ei^^’s 
and Ei^i's are the single particle energy eigenvalues ob- 
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tained by diagonalizing the up and down spin Hamilto¬ 
nians and respectively. 


B. Energy band structure 

To make this present communication a self contained 
study let us first start with the energy band structure of 
a finite width zigzag nano-ribbon. 

Non-interacting case: To establish the energy dispersion 
relation of a zigzag nano-ribbon we find an effective dif¬ 
ference equation analogous to the case of an infinite one¬ 
dimensional chain. This can be done by proper choice of 



FIG. 4: (Color online). Unit cell configuration of a zigzag 
nano-ribbon. 

a unit cell from the nano-ribbon. The schematic view of 
the unit cell configuration with Ny pairs of B-A atoms 
in a zigzag nano-ribbon is shown in Fig. With this 
arrangement, we express the effective difference equation 
of the nano-ribbon in the form, 

{El - Sa)'ipj^a = Ttpj-^l,a -b 7^V'i-l.cr (5) 


where. 


/ 




' 4 ^ jlB^a 
4^jlA,(T 
4’j2B,a 


\ 


V 4^jNyA^(7 / 


( 6 ) 


£ and T are the site-energy and nearest-neighbor hopping 
matrices of the unit cell, respectively. I is a {“^Ny x 
2iVy) identity matrix. According to our convention, the 
translational invariance of the nano-ribbon exists along 
the x-direction and we can write in terms of the 
Bloch waves and then Eq. takes the form, 

{El - £„) = (7) 


where, A = ^/3« is the horizontal separation between two 
filled magenta or yellow circles situated at two successive 


unit cells, a is the length of each side of a hexagonal ben¬ 
zene like ring. Finally, we solve Eq. [^to get the desired 
energy dispersion relation {E vs. kx) of the ribbon. 

As illustrative example in Fig. we display the varia¬ 
tion of energy levels (blue curves) as a function of wave 
vector kx for a finite width zigzag nano-ribbon consider¬ 
ing Ny = 4. At A = 0, nearly flat bands appear in the 
spectrum. The electronic states corresponding to those 
almost flat bands are strongly localized near the zigzag 
edges of the tube. The existence of these edge st ates hav e 
also been reported earlier by some other groupj^l ^^l^^ l 

With this energy band structure of a finite width nano¬ 
ribbon we now pay attention on the variation of energy 
levels of a nanotube. In order to get the nanotube from 
the ribbon we apply periodic boundary condition along 
the x-directiorP^ which results quantized values of kx, 
and the quantized wave numbers are expressed from the 



FIG. 5: (Color online). Energy levels (blue curve) as function 
of kx for a finite width zigzag nano-ribbon considering Ny = 4. 
The discrete eigenvalues (filled black circles) of a nanotube 
with Nx = 12 and Ny = 4, in absence of AB flux cj), are 
superimposed. Here we set U — 0. 


relation kx = ‘iirUx/NxE, where nx is an integer lies 
within the range: — A(j;/4 < Ux < A(r/4. Plugging the 
quantized values of kx in Eq. we can easily determine 
the eigenvalues of a finite sized nanotube. In Fig. we 
represent the discrete energy eigenvalues (filled black cir¬ 
cles) for a zigzag nanotube considering Nx = 12 and 
Ny = 4. It is to be noted that the results displayed in 
Fig. [^correspond to the case when AB flux (j) is set equal 
to zero. With these parameter values of the nanotube 
kx gets six quantized values (—tt/A, —27r/3A, —tt/SA, 0, 
7r/3A and 27r/3A), and therefore, total 48 energy values 
are obtained since Ny is fixed at 4. 

Interacting case: In the presence of e-e interaction energy 
levels get modified significantly depending on the filling 
of electrons. The results calculated for a particular value 
of U are presented in Fig. [^ where we set Ny = 4. In 
the half-filled band case, a gap opens up at the Fermi 
energjP^ which is consistent with the DFT calculation^^ 
and the gap increases with the value of U. A careful 
investigation also predicts that the full energy band gets 
shifted by the factor U/2. 
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FIG. 6: (Color online). Energy levels as function of for 
a finite width zigzag nano-ribbon considering Ny = 4 and 
U = 1.4, where (a) and (b) correspond to the one-third- and 
half-filled cases, respectively. 

C. Energy-flux characteristics 

In this sub-section we examine the energy-flux charac¬ 
teristics of a zigzag nanotube. The results for a zigzag 
nanotube considering = 10 and Ny = 7 are dis¬ 
played in Fig. where (a) and (b) represent U = 0 
and U = 1.5, respectively. In the absence of Hubbard 



FIG. 7: (Color online). Energy-flux characteristics of a half- 
filled zigzag nanotube with = 10 and Ny = 7. (a) [/ = 0 
and (b) U = 1.5. 

interaction {U = 0), energy levels are obtained simply by 
diagonalizing the non-interacting Hamiltonian and the 
nature of the energy spectrum becomes independent of 
the total number of electrons Ne in the system. While, 


for the interacting case {U ^ 0) we employ a mean-field 
scheme where the interacting Hamiltonian (Eq. are de¬ 
coupled (for a particular filling) into two non-interacting 
Hamiltonians corresponding to up and down spin elec¬ 
trons and then diagonalize both and to get the 
energy eigenvalues of the system. Since in our case we 
set Nx = 10 and Ny = 7, we get total 70 independent 
energy levels and due to their overlaps individual energy 
levels are not clearly distinguished from the spectra given 
in Fig. For identical filling factor of up and down spin 
electrons the energy levels are exactly similar both for 
and (see Fig. [^b)), and therefore, one energy 
spectrum cannot be separated from the other. At E = 0, 
the energy levels become almost flat for a wide range of 
(/), and, near (j) = ±(^o/2 they vary slowly with (j) as shown 
in Fig. [^a). In Fig. [^b) the variation of energy levels 
with (j) for a zigzag nanotube with the same parameter 
values sated above are plotted considering Hubbard in¬ 
teraction. Here we choose U = 1.5. Both for the up 
and down spin Hamiltonians the eigenvalues are exactly 
identical and they overlap with each other. In absence of 



FIG. 8: (Golor online). Energy gap (AE) as a function of 
on-site Hubbard interaction strength U for a zigzag nanotube 
with Nx = 10 and Ny = 7 in the half-filled band case when (j) 
is set at 00 / 2 . 

interaction there is no gap in the center (around E = 0) 
of the band spectrum (Fig. [^a)). As soon as the in¬ 
teraction is taken into consideration a finite gap opens 
up in the band center (Fig. [Wb)) and it increases with 
the strength of interaction (U) as shown in Fig[^ This 
energy gap is consistent with the energy gap obtained 
in the E-kx diagram (Fig. |^b)). All these energy levels 
vary periodically with 0 providing 0o (= 1 in our chosen 
unit) flux-quantum periodicity. At half-integer or integer 
multiples of (po, energy levels have either a maximum or 
a minimum (see Fig. [^. Both the energy spectra take a 
complicated look as there are many crossings among the 
energy levels particularly in the regions away from E = 0 

(Fig.g. 

In Fig. we present the variation of ground state en¬ 
ergy of a carbon nanotube with zigzag edges as a func¬ 
tion of magnetic flux 0 for the half-filled band case. Here 
(a)-(d) represent the four different cases corresponding 
to four different values of electronic correlation strength 
U = 0, 0.5, 1 and 1.5, respectively. The energy levels 
evince one flux-quantum periodicity, as expected, and 
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their energies get increased with U. 


For a particular eigenstate \ipp,a) persistent current in 
n-th channel becomes, 


D. Evaluation of persistent current in the second 
quantized form 


Following the second quantized formulatiorP^IHSl^ es¬ 
timate persistent current in individual zigzag paths of a 
nanotube, threaded by an AB flux 4>- This is an ele¬ 
gant and nice way of studying the response in separate 
branches of any quantum network. 

At first, we express the basic equation of current oper¬ 
ator Jct corresponding to spin a in terms of the velocity 
operator v„ {=x,^) as, 

I<T = 


where, x^r is the displacement operator. The velocity 
operator is obtained from the relation, 

'^fT - TT [^CTl . (9) 

ih 

We use this expression to find the velocity operator of an 
electron with spin cr in a zigzag channel n (say) in the 
form. 






-ie 


^ln,n,(T^rn+l,n,rT^ 


P 

n,£T 



\ '^n,cr | 


( 11 ) 


where, the eigenstate \'4’p,a) is written as. 


m,n 

+P?n+I,njrn + 1, n, a) + u + 1, cr)) , 

( 12 ) 


where, \m, n, a)’s are the Wannier states and ^ and 
Pmna^ are the corresponding coefficients. Simplifying 
Eq. we get the final relation of persistent charge cur¬ 
rent for n-th zigzag channel as, 


P = 


iet 

hN^ ^ 

^ m 
* RP 


( RP * fyP P 
\^m+l ,n,cr^m,n,cr^ 


-iO 


-r\P * -B'^ ry^ p 

^m,n,crA-^m+l,n,£T^ A^m—1 ,n,(T^m,n,cr^ 

J^ryP * RP p-rSf 

' ^m,n,crA^m—l,n,cr^ ) ' 


(13) 


With the above prescription we can also evaluate persis¬ 
tent current in individual armchair paths (along y direc¬ 
tion) of the nanotube. Finally, it takes the form, 




it 

2hN„ 






E (fyP * RP - RP * rrryP p*®! -I- (RP * P p“*® - P * 

V m,n,crA^m—l,n,cr^ A^m—l,n’^ ^m,n,(7^ ) ' / ^ \r'm,n,(7^m—l,n,a^ l,n,cr 


n=l,3,. 


n=2,4,. 

AAp-1 


RP _L * B^ ^ _L (RP* r,P - * RP \ 

l^m,n,<7^ ' f^m—l,n,(7^m—l,n—l,(T l,n—l,crA^m—l,n,(7/ ' / ^ V^7n,n+1 ,cr^m,n,(7 ^m,n,crA^m,n+l,(7/ 

(14) 


2=2,4,, 


where, an armchair channel (m— 1 ,to) is constructed by 
(to — l)-th and m-th lines according to our indexing. It is 
noteworthy to mention that all the calculations are done 
at absolute zero temperature (T = OAT). Now, the net 
persistent current driven by electrons of spin cr in a par¬ 
ticular channel n for a nanotube described with Fermi 
energy Ep can be determined by taking the sum of indi¬ 
vidual contributions from the lowest energy eigenstates 
upto the Fermi level. Therefore, we have, 

Inp = J2<<^- ( 15 ) 

P 


Taking the contributions from all possible channels n, 
both up and down spin (cr) electrons, the total persistent 
current in the nanotube can be expressed as, 

At = 1^4,.. (16) 

n,(T 

To judge the accuracy of the persistent current calculated 
from the present scheme we determine persistent current 
in some other ways as available in literature. Proba¬ 
bly the simplest way of determining persistent current is 
the case where first order derivative of ground sta te en- 
ergy with respect to AB flux cj) is taken into accounlP^*^. 












Therefore, we can write, 




(17) 


where, Eq{4>) is the total ground state energy for a par¬ 
ticular electron filling. But, in our present method, the 
so-called second quantized approach, there are some ad- 



FIG. 9: (Color online). Ground state energy level as a func¬ 
tion of 0 for a zigzag nanotube in the half-filled band case 
considering = 20 and Ny = 8. (a), (b), (c) and (d) corre¬ 
spond to 17 = 0, 0.5, 1 and 1.5, respectively. 

vantages compared to other available procedures. With 
the second quantized formulation we can easily determine 
current in any branch of a complicated network and the 
evaluation of individual responses in separate branches 
helps us to elucidate the actual mechanism of electron 
transport in a more transparent way. 


E. Current-flux characteristics 

Now in this subsection we focus our attention on the 
behavior of persistent current in a zigzag nanotube. Here, 
also we adopt the unit where c = /i = e = l, fixt = —1 
and measure all the physical quantities in unit of t. 

To illustrate the behavior of persistent current in sep¬ 
arate branches of a zigzag carbon nanotube we show in 
Fig. the variation of persistent current in individual 
zigzag paths as a function of flux 0 for the half-filled 
case considering = 20 and Ny = 7, where (a)-(g) 
correspond to lst-7th zigzag channels of the tube, re¬ 
spectively. The Hubbard correlation strength U is fixed 
at 1.2. Currents carried by up and down spin electrons 
are displayed in each of these figures simultaneously and 
they are exactly superposed on each other as the mag¬ 
nitudes and behaviors are exactly similar to each other. 
All these current profiles exhibit 0o flux-quantum peri¬ 
odicity. There are no similarity in magnitudes of per¬ 
sistent currents corresponding to separate channels but 
also another similarity carried by them. Interestingly we 


observe that is exactly identical to and, similarly 
for the (/ 2 ^, Iq^) and (/a^, pairs. the current 

in the middle channel, becomes the isolated one since we 
have chosen Ny = 7. This is true for any zigzag nan¬ 
otube with odd Ny. For a tube with even Ny, currents 
are pairwise identical. Summing up the individual cur¬ 
rents in seven zigzag channels we get separately the net 
persistent current carried by up and down spin electrons 
in the nanotube which is presented in Fig. [Mh) and 
the total current is obtained by taking both contribu¬ 
tions Irt and Iti and it is depicted in Fig. [To|(i). Now 
the total current derived from the conventional method 
where first order derivative of the ground state energy 
is taken into account, is displayed in Fig. [m These two 
net currents calculated from the two different schemes are 
exactly identical to each other. It emphasizes that the 
net contribution of persistent current in a zigzag carbon 
nanotube comes only from the individual zigzag channels, 
not from the armchair paths. To justify it in Fig. |12| we 
present the variation of persistent current in an armchair 
path as a function of 0 for a half-filled zigzag nanotube 
considering N^ = 20 and Ny = 7, which clearly shows 
zero current for the entire range of 0. 

In Fig. we have already observed a few almost flat 
energy levels. These energy levels contribute a very little 
to the persistent current while, on the other hand, the 
energy levels with larger slopes provide large persistent 
current. Also, at the minima or maxima points persis¬ 
tent current becomes zero which is quite obvious since the 
current is obtained by taking the first order derivative of 
the eigen energy with respect to flux 0 (Eq. . This pe¬ 
culiar nature of the energy levels invokes the current am¬ 
plitude to become filling dependent and we elaborate it in 
this section. To clarify this feature in Fig.j^the current- 
flux characteristics of a zigzag nanotube, with N^ = 14, 
Ny = 6 and U = 1, is depicted where four different fig¬ 
ures correspond to the four different cases of electron 
fillings, Nf, = 10, 20, 30, and 82. In all these cases, per¬ 
sistent current varies periodically with flux 0, exhibiting 
00 flux-quantum periodicity. Also, we observe that there 
are multiple kinks in the current profiles at different val¬ 
ues of 0 when the number of electrons remains lower than 
that required for half-filling of the band. These kinks are 
associated with the multiple crossings of energy levels, as 
shown in Figs. |13[ a)-(c) where the number of electrons 
are respectively, 10, 20, and 30. This is quite analogous 
to the feature of persistent current observed in conven¬ 
tional multi-channel mesoscopic cylinders. The behavior 
of persistent current gets significantly modified when the 
nanotube becomes half-filled or nearly half-filled. As il¬ 
lustrative example Fig. HlKd) is depicted. As we have 
already shown a few results for half-filled band case, here 
we choose Ng = 82 i.e., the nanotube is very near to the 
half-filled band condition. For this choice of parameter 
values it is examined that all the kinks disappear making 
the variation of persistent current quite smoother. This is 
analogous to the behavior of persistent current observed 
in traditional single-channel mesoscopic rings. For the 
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FIG. 10: (Color online). Persistent current in individual zigzag paths as a function of (j) for a half-filled zigzag nanotube 
{Nx = 20 and Ny = 7) with U = 1.2, where, (a)-(g) correspond to lst-7th zigzag channels of the tube, respectively. The net 
current corresponding to both up and down spin electrons are displayed in (h) while in (i) total persistent current is shown. 



FIG. 11: (Color online). Total persist ent current obtained in 
a traditional derivative approach (Eq. |17[ ) as a function of (f) 
for the same parameter values mentioned in Fig. |10| 


FIG. 12: (Color online). Persistent current in an armchair 
path as a function of cji for ^ zigzag nanotube with the same 
parameter values as mentioned in Fig. |10| 


cases when the nanotube is far away from half-filling, 
current amplitudes are quite comparable to each other 
(see Figs. [^a) -(c)). On the other hand, when the tube 
is nearly half-filled current amplitude remarkably gets 
suppressed and this suppression is very much clear from 
Fig. ^d). The reason behind this enormous reduction 
of current amplitude can be understood clearly when we 


look at the characteristics of Fig. At half-filling 
or very close to half-filling, the top most filled energy 
level lies in the nearly flat region i.e., around E = 0 (see 
Fig.j^a)) and it contributes a little to the current. More¬ 
over, when [/ 0, there is gap in the midband region. 

Now, for a particular filling we find the net persistent 
current taking the sum of individual contributions from 














10 



FIG. 13: (Color online). Current-liux characteristics of a 
zigzag nanotube with Nx = 14 and Ny = 6, where (a), (b), (c) 
and (d) correspond to Ns = 10, 20, 30 and 82, respectively. 
U is fixed at 1. 


the lowest filled energy levels, and, in this process only 
the contribution which comes from the highest occupied 
energy level finally survives and the rest disappear due to 
their mutual cancellations. It leads to the enormous re¬ 
duction of persistent current amplitude in the half-filled 
or nearly half-filled band case. This feature is indepen¬ 
dent of the size of the nanotube. 

Before we end this section we can emphasize that, the 
current amplitude in a zigzag nanotube is highly sensitive 
to the electron filling and this phenomenon can be uti¬ 
lized in designing a high conducting to a low conducting 
switching operation and vice versa. 


III. AN ORDERED-DISORDERED SEPARATED 
NANOTUBE IN PRESENCE OF A MAGNETIC 
FLUX 


A. The model 

The schematic diagram of the model quantum ordered- 
disordered separated nanotube is illustrated in Fig. [T^ 
where M co-axial rings are vertically attached and each 
ring contains N atomic sites. The cylinder is subjected 
to an AB flux (j). For a comparative study we will explore 
the behavior of persistent currents in a fully disordered 
and an ordered-disordered separated cylinder simultane- 
ouslji^. In order to get an ordered-disordered separated 



FIG. 14: (Golor online). Schematic diagram of a mesoscopic 
cylinder having 6 rings where each ring contains 8 atomic 
sites. 

cylinder disorders are introduced into M/2 number of 
rings of the lower half of the system (green sites) whereas 
the remaining part (magenta sites) becomes the ordered 
one. On the other hand, full disordered system is ob¬ 
tained by putting random disorder in all the atomic sites. 
A non-interacting single-band tight-binding (TB) Hamil¬ 
tonian is used to describe the model quantum system 
which reads as. 


This section illustrates the behavior of persistent cur¬ 
rent in an ordered-disordered separated mesoscopic cylin¬ 
der. Persistent current in conventional disordered sys¬ 
tems have already been reported by many physicists and 
Anderson type localization^ due to disorder is not a new 
one. But, recent breakthroughs in the growth of semicon¬ 
ductor nanowires have opened up great opportunities to 
revolutionize technologies in nanoscale electronic^^. Re¬ 
cent experiments on nanowires have already yielded some 
results in contrast to the phenomenon of strong localiza¬ 
tion due to doping. For instance, Cui et alW^ reported 
that the carrier mobility in boron-doped and phoporus- 
doped silicon nanowires under low dopant concentration 
is extremely low compared to bulk silicon but the con¬ 
ductance becomes diffusive and about five orders of mag¬ 
nitude larger with heavy doping. This motivates us to 
investigate whether these kinds of shell doped nanowires 
give rise to new features in the context of persistent cur¬ 
rent or not. 


^ ^ ^ T ty{e “b h.C.^ 

m,n 

+ty{e"^'‘cl^^^Cm+i,n + h.c.)] (18) 

where, Cm.n is the on-site energy. For the impurity sites, 
the site energies Cm.n are selected randomly from a ‘Box’ 
distribution function of width W = 1. ty and refer to 
the intrachannel and interchannel nearest-neighbor cou¬ 
plings, respectively. Here, {m,n) is the co-ordinate of a 
lattice point where, m and n run from 1 to M and iV, re¬ 
spectively. To incorporate the effect of magnetic flux we 
consider 6r = 2Tr(j)/N with the threading magnetic flux 
(j) measured in unit of the elementary flux-quantum c/q. 
The phase factor along the vertical direction is taken 
as zero, cl^ „ {cm,n) is the creation (annihilation) oper¬ 
ator of an electron at the site {m,n). At absolute zero 
temperature, the total persistent current of the cylinder 
is determined from the relation of Eq. [m where Eq is 
the ground state energy of the system which is obtained 
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by taking the sum over lowest N^. (number of electrons) energy levels are much more different than that in the 
energy levels. case of full disordered system. In order-disordered sepa¬ 

rated nanotube the energy levels become flatter mainly 
from the boundary of the spectra. 

B. Band structure 


In this sub-section we briefly discuss the energy band 
spectra of a mesoscopic cylinder. To find an analytic form 
of the energy levels we take a very small sized cylinder 
having only two rings with 10 sites in each ring. For this 
smallest possible size of the cylinder, expression of the 
energy levels takes the form, 


E - Ety -j- ‘ItyCOS 


27r(n -I- (j)) 

Ns4>o 


(19) 


where, iVg(= MxN) is the total number of atomic sites 
of the cylinder. In Fig. we plot the variation of energy 
levels as a function of magnetic flux (j). The numerical 
values of the parameters are taken as ty = ty = —1 and 
a = 1, where a is the lattice parameter. Two energy 
bands consisting of N number of energy levels within the 
boundary E' = —1 to -1-3 and —3 to -1-1 are found to 
overlap in the region E = —1 to -1-1 (in unit of ty). The 
width of the overlap region can be tuned by tuning the 
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FIG. 15: (Color online). Energy-flux characteristics for a 
mesoscopic cylinder with 2 rings where each ring contains 10 
atomic sites. 


vertical hopping strength {ty). Each energy level exhibits 
one flux-quantum quantum periodicity. The variation of 
the energy levels of an ordered multi-channel system are 
already discussed in the previous section. Effect of the 
disorder can be attributed to the removal of degeneracy 
of the energy levels which can be realized from Fig. |16| 
The results are calculated numerically for a particular 
disorder configuration. In presence of impurity gap opens 
up at the crossing points and width of this gap increases 
with the increase of the disorder strength making the 
slope of the energy levels smaller and smaller. This makes 
a remarkable change in magnitude of persistent current. 
We discuss it in the forthcoming sub-section. 

Not only the full disordered cylinder here we also 
present the energy spectra of an ordered-disordered sep¬ 
arated cylinder in Fig. 17 where (a) and (b) correspond 
to two different values of the disorder strengths, W = 2.5 
and 5, respectively. Here, the nature of the variation of 


C. Results and discussion 

In order to study the effect of disorder on persistent 
current in Fig. we plot the typical persistent current 
amplitude / at a particular value of <() as a function of the 
disorder strength W for an ordered-disordered separated 
cylinder. Figure (a) corresponds to the results for a half- 
filled cylinder consisting of 10 rings with 12 sites in each 
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FIG. 16: (Color online). Energy-flux characteristics for a 
mesoscopic cylinder with the same parameter values men¬ 
tioned in Fig. in presence of disorder at all atomic sites 
{W = 3.5). 

ring and the magnetic flux (j) = 0.4. On the other hand, 
figure (b) represents the results for a half-filled cylinder 
having 8 rings (10 atomic sites in each ring) in presence of 
magnetic flux (j) = 0.2. Here we compute the root mean 
square of the current amplitude taking the average over 
30 random disordered configurations. We observe that 
for a full disordered cylinder the current amplitude de¬ 
creases with the rise of impurity strength W (red curve), 
while in the case of ordered-disordered separated cylinder 
the current amplitude initially decreases upto a certain 
value of W after which it slowly increases with the rise of 



FIG. 17: (Golor online). E-d> curves for an ordered-disordered 
mesoscopic cylinder with 2 rings and 10 atomic sites in each 
ring. 

the impurity strength (green curve). This phenomenon 
can be explained as follows. In the presence of disorder 
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current amplitude decreases due to the localization of the 
energy levels which is the so-called Anderson localization. 
The more impurity strength results the more reduction 
in current amplitude. This is true only for a full disor¬ 
dered system. But, for the ordered-disordered separated 
system there is a different physical picture. As we gradu¬ 
ally tune the disorder strength towards the higher value, 
the ordered and disordered regions get decoupled from 
each other and after a critical limit of W, (say Wc), the 
cylinder behaves as composed of two completely decou¬ 
pled regions. At this situation only the ordered region 
contributes to the current, and thus, the current ampli¬ 
tude gets enhanced with the disorder strength W. The 
critical value Wc depends on the system size and other 




w 


w 


FIG. 18: (Color online). I-W characteristics of an ordered- 
disordered separated cylinder in the half-filled band case, (a) 
M = 10, A = 12 and (j> = 0.4 and (b) M = S, N = 10, and 
0 = 0.2. The red and green color correspond to full disordered 
and half-disordered cylinder. 


parameter values like magnetic flux. As for example, it 
is 3 and 2.5 for two different system sizes as shown in 
Figs. [T8|(a) and (b), respectively. 

To summarize, in this section we have explored the 
behavior of persistent current in a disordered and an 
ordered-disordered separated cylinder in presence a mag¬ 
netic flux. Most interestingly, we see that the current 
amplitude shows an anomalous behavior with the in¬ 
crease of impurity strength for an ordered-disordered sep¬ 
arated cylinder in contrast to the completely disordered 
one. This study may be helpful for exploring localization- 
delocalization transition in shell-doped nanotubes. 


IV. A BINARY ALLOY RING WITHOUT 
EXTERNAL ELECTRODES 


In this section we undertake an analysis of the band 
structure and persistent current in an isolated binary al¬ 
loy ring (no source and drain electrodes) enclosing a mag¬ 
netic flux ((|S3 


A. The model 

Let us concentrate on the simplest model of a binary 
alloy ring as shown in Fig. [T^ where the ring consists 
of two different types of atoms placed alternately in a 
regular pattern. They are characterized by two different 
on-site potential energies, namely, a and /?. The ring 
is subjected to an AB flux (j). Within a non-interacting 
single-band TB framework we illustrate the model of bi- 


tr 
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FIG. 19: (Golor online). A binary alloy ring, threaded by a 
magnetic flux (j), is composed of two different types of atomic 
sites, viz, a and {3 those are represented by filled blue and red 
circles, respectively. 

nary alloy ring and the TB Hamiltonian reads, 

Hr = ^{eicjci ( 20 ) 

i 

where, the on-site energy, ei takes two values Cq, and 
corresponding to two different sites a and /?, respectively. 
tr is the nearest-neighbor hopping integral. The phase 
factor 9 = 2TTcj)/N of the Hamiltonian takes an account 
of the effect of the magnetic flux </) threaded by the ring 
which is measured in unit of the elementary flux-quantum 
(j)Q. cj (ci) is the creation (annihilation) operator of an 
electron at the site 1. Here, I runs from 1 to N, where N 
is the total number of sites in the binary ring. 


B. Energy spectrum 


Before addressing the main points i.e., the character¬ 
istic features of persistent current in an ordered binary 
alloy ring, let us have an idea about the energy band 
structure of the system. The analytical expression of en¬ 
ergy dispersion relation for the ordered binary alloy ring 
is as follows. 


E = 


Ect + 6/3 
2 



2 

-b 4^2 cos'^{ka) 


( 21 ) 


where, a is the lattice spacing and k is the wave vec¬ 
tor. The periodic boundary condition of the ring sets the 
quantized values of k in presence of the AB flux (j) as, 


k = 


2tt 

Ao 



( 22 ) 
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where, n is an integer and it is restricted within the range: 
—N/2 < n < N/2. Throughout our manuscript we con¬ 


sider Ca = —ep = e and Eq. 21 is modified according to 
this condition as, 


E = ±-\/ -I- 4^2 cos'^{ka). 


(23) 


In Fig. 20 we plot the energy levels as a function of flux 
(p, obtained from Eq. for a 40-site binary alloy ring 
considering e = 1 and = 1. Looking at Fig. two 

different sets of energy levels are noticed to form two 
quasi-bands separated by a finite energy gap. This gap, 
on the other hand, is tunable by the parameter values de¬ 
scribing the TB Hamiltonian Eq. The origin of two 
different sets of energy levels is also clearly understood 


from Eq. ^ The energy levels of Fig. have either a 
maximum or a minimum at half-integer or integer mul¬ 
tiples of flux-quantum and it results vanishing nature of 



FIG. 20: (Color online). Energy-flux characteristics of an 
ordered binary alloy ring {N = 40) considering Cc = —ep = 1 
and tr = 1. 00 is set at 1. 

persistent current at these specific values of 0, since the 
current is obtained by taking the first order derivative 
of energy E^p) with respect to flux 0. All these energy 
levels vary periodically providing 0o flux-quantum peri¬ 
odicity, 00 being 1 in our chosen unit (c = /i = e = 1). 


C. Persistent current 

Our task of calculating persistent current for individual 
energy eigenstates is now easier as we know the energy 
eigenvalues of the ring as a function of flux 0. It is simply 
the first order derivative of energy with respect to flux 
0. Therefore, for an n-th energy eigenstate we can write 
the expression for the current as, 

I^ = ±( sin[^^{n + <P/(t>o)] 

where, -l-ve or —ve sign appears in the current expression 
depending on the choice of n i.e., in which sub-band the 
energy level exists. To get total persistent current I for a 


particular filling Afg, we take the sum of individual con¬ 
tributions from the lowest Ng energy eigenstates as we 
do our calculations at absolute zero temperature. The 
expression is, 

AT, 

I=Y.In (25) 

n—1 


Following this way, we plot the variation of persistent cur¬ 


rent for an ordered 120-site binary alloy ring in Fig. 21 


The values of the parameters considered in this figure are 
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FIG. 21: (Golor online). Current-flux characteristics of an 
ordered binary alloy ring {N = 120) considering = —ep = 1 
and tr = 1. The blue and purple colors in (a) correspond to 
Nr = 15 and 35, while in (b) they correspond to Ne = 10 
and 30, respectively. The lattice spacing a is set at 1 and we 
choose 00 = 1. 


e = 1 and tr = 1. In both figures, Fig.j^a) and (b), the 
current profiles show saw-tooth like variation as a func¬ 
tion of flux 0, similar to that of traditional single-channel 
mesoscopic ring^. The purple and blue colors corre¬ 
spond to two different fillings of the band. For odd num¬ 
ber of electrons the results are shown in Fig.[2l|(a) while, 
in Fig. [13b) the results are given for even number of 
electrons. The sharp transitions at half-integer (for odd 
Ng) or integer (for even Ng) multiples of flux-quantum 
(0o) in persistent current appears due to the crossing of 
energy levels at these respective values of 0. Quite inter¬ 
estingly, we also examine that the current shows always 
diamagnetic response irrespective of the filling factor. 


V. A BINARY ALLOY RING WITH EXTERNAL 
ELECTRODES 

Upto now we have discussed the basic features of per¬ 
sistent current in different isolated mesoscopic systems. 
In this section we extend our analysis of persistent cur¬ 
rent to an open system where a binary alloy rinj^ is 
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clamped between two semi-infinite one-dimensional (ID) 
electrodes. The behavior of persistent current in presence 
of transport current will also be analyzed. 


A. The model 


Let us start by referring to Fig. where a binary 
alloy ring, threaded by a magnetic flux is attached 
to two semi-infinite one-dimensional metallic electrodes, 
namely, left-lead and right-lead, via two atomic sites la¬ 
beled as /i and v. The total numbers of identical pairs of 
sites in the binary ring is Ni. Now we incorporate 
some foreign atoms denoted by 7 having on-site poten¬ 
tials e-y in any one of the two arms of the ring. Here, we 
choose the upper arm in this purpose. There are N 2 iden¬ 
tical foreign atomic sites, embedded together in a small 
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FIG. 22: (Color online). A binary alloy ring in presence of 
some identical foreign atoms (labeled as 7 sites), threaded by 
an AB flux cj), is attached to two semi-infinite one-dimensional 
metallic leads. ^ and v are the connecting sites. 


portion of the ring. These 7 sites are often referred to 
as impurity sites in the present manuscript. We enumer¬ 
ate the atomic sites of the two side-attached leads in a 
particular way, as shown in Fig. |22| A single-band non¬ 
interacting TB framework is used to describe the entire 
system. For the full system we can partition the total 
Hamiltonian as a sum of three terms like, 

H = Hr + Hl + Ht (26) 

where, Hu, and Ht represent the Hamiltonians for 
the ring, leads (left and right) and coupling between the 
ring and leads, respectively. The ring Hamiltonian {Hji) 
takes the form exactly similar to Eq. but instead of 
two possible on-site potentials here e; has three possibil¬ 
ities for three different atomic sites (a, /3 and 7 ). By 
means of some external gate voltage Vg, the site energy 
e-y can be tuned and accordingly the site-energies are 
changed. Thus we express the on-site energy of a sin¬ 
gle 7 atom like, = e^ + Vg, where e° is the site energy 
in absence of any external potential. The other two terms 
of Eq. |26[ H^ and Ht, can also be written in a similar 


fashion as, 

Hl = io + h.C.^ +tQ (^bl^bm+1 + h.c. 


m<0 


m>l 


left lead 


right lead 


and, 


(27) 


Ht = 


(^Lblcg, + Tiib\c„^ + h.c. (28) 


where bl^ (bm) are the creation (annihilation) operator 
of an electron at the site m of the leads and to represents 
the nearest-neighbor hopping strength within these leads. 
We set the site energy for the identical sites in the leads 
to zero and due to this reason in Eq. [^we have omitted 
the site-energy terms. Here, tt is the coupling strength 
between the left lead and the ring, while it is for the 
other case. 


B. Wave-guide theory 

To find transmission probability across the ring and 
also to calculate persistent current in such an o pen ring 
geometry we adopt the wave-guide theorjS^EH in the 
present sub-section we describe the formulation very 
briefly. 

We begin with the Scrodinger equation i7|'(/') = E\ip), 
where {ip) is the stationary wave function of the entire 
system. In the Wannier basis it {\tp)) can be expressed 
as, 

\pj)=Y,Bm\m)+Y,Bm\m) + Y,Ci\l) (29) 

m<0 m>l I 

left lead right lead ring 

where, the co-efficients Bm and C; correspond to the 
probability amplitudes in the respective sites. Keeping 
in mind the periodicity of the ordered binary ring we 
write the wave functions associated with the electrons as 
a plane wave and the wave amplitudes in the left and 
right leads are, 

= +re-^^^, for to<0 (30) 

and 

Bm = for m > 1 (31) 

where, r and t are the reflection and transmission ampli¬ 
tudes, respectively, k is the wave number and it is related 
to the energy E of the incident electron by the expression 
E = 2to cosk. The lattice spacing a is set equal to 1. 

In order to find out the transmission amplitude t, we 
have to solve the following set of coupled linear equations. 

EBo = toB-i tlC^, 

{E-ei)Ci = UE^Ci+i+tre-^’^Ci-i + TLBoSi^f, 

+ trBiSi^^ 

EBi = toB2 ttiCu 


( 32 ) 
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where, the co-efhcients -Bq, B_i, Bi and B 2 can be easily 
expressed in terms of r and t by using Eqs. and 
and they are in the form: 

Bo = l + r 
B_i = — 2 isinfc 

Bi = 

B 2 = (33) 

Thus, for a particular value of E we can easily solve the 
set of linear equations and find the value of t. Finally, 
the transmission probability across the ring becomes 

T{E) = (34) 

Now, to compute the persistent current between any two 
neighboring sites in the binary ring we use the following 
relation. 


/u+i = ^Im(c/Q+ie-*^). (35) 


C. Transmission and average density of states 


Throughout our calculations we set e = 1, e° = 0, 
tr- = 1, eo = 0 and to = 2. The energy scale is measured 
in unit of tr- 

Two-terminal transmission probability T (orange 
color) as a function of injecting electron energy E for 
some typical binary alloy rings considering different num¬ 
ber of impurity sites is displayed in Fig. 23 where (a), (b) 
and (c) correspond to three different numbers of impurity 
atoms. Figure (a) represents the transmission spectrum 
of the binary alloy ring with 28 pairs of a-P atoms but 
without any foreign impurity atoms while (b) and (c) 
correspond to N 2 = 12 and 22 , respectively, keeping Ni 
fixed to 28. The average density of states (ADOS) is also 
superimposed in each spectrum. In all these cases the 
magnetic flux p is set equal to zero. The results are quite 
interesting. In absence of any impurity site the transmis¬ 
sion spectrum is characterized by two bands separated 
by a finite gap and the gap between these two bands are 
also tunable by the parameter values describing the sys¬ 
tem. Looking back at Fig. the band splitting for the 
ordered binary alloy ring is easily understood and we pre¬ 
dict that the spectrum is just the finger print of the sys¬ 
tem energy levels. This transmission spectrum exactly 
overlaps with the ADOS profile {p-E spectrum) which 
ensures that electronic transmission takes place through 
all the energy eigenstates of the binary alloy ring and 
they are extended in nature. The situation becomes re¬ 
ally interesting when some additional impurities (7 sites) 
are introduced in any part of the binary alloy ring. Due 
to the inclusion of such atomic sites some energy levels 
appear within the band of extended regions those are no 
longer extended, but they are almost quasi-localized and 
do not contribute to the electronic transmission. This 


behavior is elaborated in Figs. [2^b) and (c), where T- 
E and p-E spectra are superimposed with each other 
for two different numbers of impurity sites (A^ 2 )- In the 
band center, several energy levels appear and these lev¬ 
els do not provide any contribution to the transmission 
of electrons. They are almost localized in nature. Mainly 
the states within the two bands are responsible for elec¬ 
tron transmissions. The number of the almost localized 
energy levels increases with the increase of the number 
of the impurity sites, and for sufficiently large number 
of impurities they form a quasi-energy band of localized 
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FIG. 23: (Color online). Transmission probability (orange 
color) and average density of states (dark-blue color) for some 
typical binary alloy rings with fixed number of a-j3 pairs 
(Ai = 28), but different values of impurity sites N 2 , where 
(a) N 2 = 0, (b) N 2 = 12 and (c) N 2 = 22. For these three 
cases we set p = 1 and choose u = 29, 35 and 40, respectively. 
Other parameters are: = 0, = rij = 1 and (f> = 0. 


states. The location of the localized energy band can be 
shifted towards the edge of extended regions simply by 
tuning the site energy of these foreign atoms, and this can 
be done by means of applying an external gate voltage 
Vg. We utilize this feature to make the binary alloy ring 
behave like an extrinsic semiconductor, either p-type or 
n-type by tuning the Fermi level to an appropriate place. 

We plot Fig. to explore the semiconductor-like be¬ 
havior of the binary ring. The variation of transmission 
probability and the ADOS are shown in this figure. Here, 
we consider a binary ring with 60 identical pairs of a-/3 
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sites and 44 number of impurity sites. Two different cases 
are exhibited in Fig. 24 'a) and (b) for two different values 
of the on-site potential energies of the 7 atoms. Figure 
(a) represents the transmission spectrum of the binary 
alloy ring when is set at 0.6 while the other figure 
(Fig.[24|^b)) corresponds to e-y = —0.6. In presence of the 
impurity atoms we get almost quasi-energy bands (ADOS 
spectra) and quite interestingly we observe that when 
is fixed at 0 . 6 , a localized energy band for a wide range 
of energy is formed along the left edge of the extended 
region (Fig. [2^a)). Now, if the Fermi level is set around 
E = —1.7, then many electrons in the localized region 
below the Fermi level can jump easily, even at much low 
temperature since the energy gap is almost zero, to the 
extended regions and can contribute to the current. As a 
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FIG. 24: (Color online). Transmission probability (orange 
color) and ADOS (dark-blue color) as a function of energy for 
a binary alloy ring (Ai = 60) with 44 impurity sites (Ni = 
44), where (a) and (b) correspond to e.y = 0.6 and —0.6, 
respectively. Other parameters are: (j) = 0, tl = tr = 1.5, 
— 1 and V = 83. 


result large number of excess electrons become available 
in the conduction band region which behave as n-type 
carriers. This movement however depends on the local¬ 
ized energy levels and also the available extended energy 
states. On the contrary, when the site-energies of the 7 
atoms are tuned to another value, say, — 0.6 an exactly 
opposite behavior is obtained. In this case the wide band 
of localized states is formed in the right edge of the ex¬ 
tended region (see Fig. 24'b)). Now, if the Fermi level 
is tuned to appropriate place, around E = 1.7 in this 
case, then the electrons from the filled extended levels 
below the Fermi level hop to the nearly empty localized 
levels, and these electrons do not contribute anything to 


the current. But, the absence of electrons are realized 
by holes in the extended regions which can carry current 
and the system behaves like a p-type carriers. Before we 
end this discussion, we can emphasize that by setting the 
Fermi level in appropriate places our model quantum sys¬ 
tem can be tailor-made to use as a p-type or an n-type 
semiconductor. 

At the end of this sub-section we would like to mention 
a few points. To establish the fact that how such a geome¬ 
try can be utilized as a p-type or an n-type semiconductor 
with appropriate choice of the Fermi level we have taken 
a particular set of parameter values. As we are doing 
a model calculation it is very easy to take some specific 
values of the parameters and use in our numerical cal¬ 
culations, but all these physical phenomena are exactly 
invariant with the change of the parameter values. Only 
the numerical values will be altered. These features are 
also exactly valid even for a non-zero value of magnetic 
flux (j). The physical picture will be much more appealing 
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FIG. 25: (Golor online). Persistent current (lu) in the upper 
arm (red and violet color) as a function of energy E for a 
binary alloy ring {Ni = 20) in the absence of impurity atoms 
{N 2 = 0), where (a) <?!> = 0 and (b) (j) = i/io/d. Other parame¬ 
ters are: tl = tr = 0.8, p = 1 and = 21. For this ordered 
ring, persistent current in the lower arm (II) is exactly iden¬ 
tical to lu- ADOS profile (gray and light-blue color) is also 
superimposed in each spectrum. 


if we consider larger rings with more impurity sites and it 
gives us the confidence to propose an exp eriment in this 
way. In a recent work Bellucci et have done a de¬ 

tailed study of magneto-transport properties in quantum 
rings considering tunnel barriers in the presence of mag¬ 
netic field and shown how metal-to-insulator transition 
takes place in such a geometry. They have also estab¬ 
lished that by controlling the strength and the positions 
of the barriers, the energy shift can be done in a tunable 
way. This is quite analogous to our present study, and 
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so, an experiment in this regard will be challenging. 


D. Persistent current in presence of transport 
current 

Now we focus on the characteristic features of persis¬ 
tent current in the binary alloy ring in presence of exter¬ 
nal bias voltage. We also study the role of magnetic field 
in this context. In Fig. |25| upper arm current, Ijji in the 
upper arm of an ordered binary alloy ring as a function 
of energy E is presented where (a) and (b) correspond 
to (j) = 0 and ^io/4, respectively. Due to conservation of 
current at the junction points, current between any two 
sites are equal to each other. We refer upper arm current 
to the bond current between any two nearest-neighbor 
atomic sites of the upper arm. The ring is symmetrically 
coupled to the side attached leads i.e., the upper and 
lower arms have identical length, and accordingly, the 
current 1^ in the lower arm becomes exactly identical to 
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FIG. 26: (Color online). Persistent current as a function of 
energy for a binary alloy ring {Ni = 16) in presence of impu¬ 
rity sites {N 2 = 12) for a finite value ot (j) {(f> = 0.3), where (a) 
and (b) correspond to the currents for the upper (pink color) 
and lower (violet color) arms respectively. Other parameters 
are: Sj = 0.5, tl = tr = 0.8, /r = 1 and v — 23. ADOS 
profile (green and light-blue color) is superimposed in each 
spectrum. 


the current obtained in the upper arm. Similarly, we de¬ 
fine the lower arm current. Here, again the band-splitting 
put its mark on the current density profile. Finite current 
is available for two wide range of energies, separated by 
a finite gap, associated with the energy levels of the ring 
those are clearly visible from the ADOS profiles (green 
and light-blue color). The quite exact superposition of 
the current profile and the AVDOS prohle invokes the 
extended natures of all these energy eigenstates in our 


mind. Additionally, all the values of the current are pos¬ 
itive i.e. they are in the same phase ensuring the fact 
that this current is completely due the transport current 
between the two leads hxed at two different chemical po¬ 
tentials as there is only one driving force, the external 
bias. This phenomenon depicted in Fig.[2^a) correspond 
to zero magnetic flux. As well as the magnetic field with 
flux density (f) is switched on, not only both the positive 
and negative values of current appear indicating opposite 
phases but also they appear in regular alternative fash¬ 
ion as shown in Fig. [2^b). Moreover, the magnitudes 
are also different in these two cases indicating the effect 
quantum interference of mesoscopic regime. Therefore, 
by measuring persistent current we can directly estimate 
the nature of the current i.e., whether it is paramagnetic 



FIG. 27: (Color online). Persistent current as a function of 
magnetic flux for a binary alloy ring (Ai = 16) in presence 
of impurity sites {N 2 = 12) with Vg = 0.5 Other parameters 
are: tl = tr — 0.8, n = I and v = 23. 


or diamagnetic in nature and also predict the character¬ 
istics of energy eigenstates, which are somewhat interest¬ 
ing in the study of electron transport. Needless to say, 
magnetic field plays an important role in this context. 

Already stated in introduction the current can persists 
in this kind of loop geometry if we set the lengths of two 
arms unequal or put any kind of static disorder in path 
of the electron or in any other way we can disturb the 
wave function or create any path difference between the 
wave functions corresponding to different arms. In our 
study we incorporate impurities in the upper arm of the 
binary ring as mentioned earlier. The results are dis¬ 
played in Fig. 26 where (a) and (b) represent the cases 
corresponding to the upper and lower arms of the ring, 
respectively. To plot the curve we consider a typical bi¬ 
nary alloy ring with A^i = 16 and N 2 = 12. With the 
inclusion of impurity sites, the symmetry between the 
two arms is broken, and therefore, the currents in the 
upper and lower arms are no longer identical to each 
other as shown from the spectra (Figs. [2^a) and (b)). 
So, there must be a circulating current within the ring 
which is the so-called persistent current in presence of 
transport current. Unlike to the ordered binary alloy 
ring (Fig.|^, here all the energy eigenstates are not ex¬ 
tended in nature. Some localized energy levels appear in 
the band of extended energy states due to the presence 
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of impurity sites in the ring. This is clearly visible from 
the spectra since for these impurity levels no current is 
available. We also check the periodicity of the persis¬ 
tent current. It is shown in Fig. Persistent current 
exhibits flux-quantum periodicity. Thus, calculating 
persistent current we can emphasize the nature of energy 
eigenstates very nicely, and this idea can be utilized to 
reveal the localization properties of energy eigenstates in 
any complicated geometry. 

VI. CONCLUDING REMARKS 

To conclude, in the present article we have made a de¬ 
tailed investigation of magneto-transport in both closed 
and open mesoscopic systems. We have started with 
the discussion of a zigzag nanotube pierced by a mag¬ 
netic flux using a generalized Hartree-Fock mean field ap¬ 
proach. Based on the tight-binding model we explore the 
effect of the Hubbard interaction on the energy levels of 
the tube. After describing different numerical results we 
have established the second quantized form to evaluate 
persistent current in individual paths of a zigzag carbon 
nanotube and based on this formulation we have also pre¬ 
sented the numerical results for the distribution of persis¬ 
tent current in different branches of the nanotube. From 
the current-flux characteristics we have emphasized that 
the current amplitude in the zigzag nanotube is highly 
sensitive to the electron filling and this phenomenon can 
be utilized in designing a high conducting to a low con¬ 
ducting switching device and vice versa. Following the 


study of nanotube we have explored the persistent cur¬ 
rent in an ordered-disordered separated cylinder pene¬ 
trated by a magnetic flux. Most interestingly, we have 
seen that the current amplitude shows an anomalous be¬ 
havior with the increase of impurity strength. This study 
may be helpful for exploring localization-delocalization 
transition in shell-doped nanotubes. To discuss per¬ 
sistent current in open system i.e., system with side- 
attached electrodes we have considered a binary alloy 
ring in presence of a magnetic flux (j), as an illustrative ex¬ 
ample. Within a single-band non-interacting TB frame¬ 
work persistent current and band structure of an isolated 
ordered binary alloy ring have been analyzed. Then, we 
have explored the magneto-transport properties of a bi¬ 
nary alloy ring in presence of external electrodes. The 
effect of impurities have also been addressed. Quite in¬ 
terestingly we have noticed that in the presence of some 
foreign atoms, those are not necessarily be random, in 
any part of the ring, some quasi-localized energy levels 
appear within the band of extended energy levels. The lo¬ 
cations of these almost localized energy levels can also be 
regulated by means of some external gate voltage. This 
leads to a possibility of using such a system as a p-type 
or an n-type semiconductor by fixing the Fermi level in 
appropriate places. 

Before we end, we would like to mention that in the 
present review we have addressed some important aspects 
of quantum transport through some low-dimensional 
model quantum systems. Several other important quan¬ 
tum phenomena in such meso- and nano-scale systems 
have also been reported in recent review^^^^^^. 
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